PeV neutrinos from intergalactic interactions of cosmic rays emitted 

by active galactic nuclei 
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The observed spectra of distant blazars are well described by secondary gamma rays produced in 
line-of-sight interactions of cosmic rays with background photons. In the absence of the cosmic-ray 
contribution, these spectra would appear surprisingly hard, but the cosmic ray interactions generate 
very high energy gamma rays relatively close to the observer, and the spectra agree with the data. 
The same interactions of cosmic rays are expected to produce a flux of neutrinos with energies 
peaked around I PeV. We show that the predicted diffuse isotropic neutrino background from many 
distant sources can explain the neutrino events recently detected by the IceCube experiment. We 
also find that the flux from any individual nearby source is insufficient to account for these events. 
The narrow spectrum around f PeV implies that a typical active galactic nucleus can accelerate 
protons to EeV energies. 
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The IceCube collaboration has detected two neutrinos 
with energies 1.1 PeV and 1.3 PeV (35% systematic er- 
ror) [l[. These neutrinos are either electron or tau neu- 
trinos. The muon analysis, currently under way, is ex- 
pected to produce additional events (probably, with a 
lower energy resolution). The narrow energy range in 
which the two neutrinos have been detected indicates 
that the source has a peak somewhere around a PeV 
energy, above the experimental threshold of 0.4 PeV and 
below the Glashow resonance that enhances detector sen- 
sitivity around 6.8 PeV 0. While a number of astrophys- 
ical processes can be responsible for high-energy neutri- 
nos, only specific types of sources can produce a peaked 
spectrum around a PeV 

Narrow spectra peaked around 1 PeV were predicted to 
arise from line-of-sight interactions of cosmic rays emit- 
ted by blazars [4( if the spectrum of protons terminates 
at Sp.max 10 9 GeV 0. There is growing evidence that 
intergalactic cascades initiated by line-of-sight interac- 
tions of cosmic rays produced by active galactic nuclei 
(AGNs) are responsible for the highest-energy gamma 
rays observed from blazars 0-E3- Blazar spectra are 
explained remarkably well with secondary photons from 
such cascades 4-6], while, in the absence of such con- 
tribution, some unusually hard intrinsic spectra [l3l - [l5| 
or hypothetical new particles [lj| have been invoked to 
explain the data. Models for hard intrinsic spectra of 7 
rays can be constructed, but the natural ease with which 
secondary photons reproduce the data makes the expla- 
nation based on cosmic rays very appealing. Further- 
more, detailed theoretical models show that relativistic 
shocks can effectively accelerate cosmic rays to energies 
~ 10 s GeV, but higher energies may require rather ex- 
ceptional conditions [I3| • Secondary gamma rays explain 
the blazar spectra as long as the magnetic fields are in 
the femtogauss range flij |. 



The mechanism predicts that some neutrinos with a 
peaked spectrum should be produced in interactions of 
protons with extragalactic background light (EBL) [4j,|5|. 
We will examine whether these neutrinos can account for 
PeV neutrino events in IceCube. 

Assuming the scenario of Refs. 0-113 

, we have consid- 
ered two possibilities for the origin of IceCube neutrinos: 
a single nearby source, and a combined contribution of 
distant sources. Two important constraints have to be 
satisfied: the gamma-ray background and the cosmic-ray 
spectrum should not exceed the observed fluxes. We do 
not assume that the cosmic ray spectrum up to ultrahigh 
energies is explained by the same sources; as was pointed 
out in Ref. [19|, such a scenario disagrees with the data. 
Furthermore, we do not consider neutrinos produced in- 
side cosmic-ray sources as in Ref. [20( | . 

The total neutrino flux (summed over flavors) is un- 
certain, and the numerical value depends on the shape of 
the spectrum. For crude estimates, one can use 
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where E v is the neutrino energy. However, this numer- 
ical value can be misleading. In comparing our model 
with the IceCube data, we do not rely on this number, 
but we calculate the expected number of events using 
the spectral shape predicted by our model, as described 
below. 

For the case of one or a few nearby point sources, we 
have not been able to find an acceptable explanation of 
the IceCube events. Indeed, the neutrino required flux 
is an order of magnitude greater than the predicted flux 
from a single source shown in Ref. 0. A single source 
at a smaller distance from Earth than the blazars con- 
sidered in Ref. would produce an unacceptably large 
flux of cosmic rays. A careful numerical analysis confirms 



2 



this conclusion. We, therefore, proceed to considering 
the second possibility: a diffuse background from distant 
sources. 

For the diffuse flux calculation we use the numerical 
code described in detail in Ref. |2lj] . The code is based 
on kinetic equations; it calculates the propagation of nu- 
cleons, stable leptons and photons using the standard 
dominant processes, i.e. pion production by nucleons, 
pair production by protons and neutron /3-decays. 
For electron-photon cascade development, it includes 
pair production and inverse Compton scattering. We also 
take into account neutrino oscillations on their way from 
the site of production to the observer. Since the distance 
traveled by neutrinos is much greater than the oscillation 
length, muon neutrinos oscillate into tau neutrinos with 
a 50% probability. The resulting spectrum has a flavor 
ratio of approximately (1:1:1). We note that in our nu- 
merical calculations we use the actual mixing angles in 
the tri-bimaximal neutrino mixing approximation. 

A number of different models have been advanced for 
EBL [iirEl. There are some upper bounds on EBL 
in the literature that were based on observations of dis- 
tant blazars but were derived without taking into ac- 
count the cosmic ray contribution. When the cosmic rays 
are included, these bounds on EBL are relaxed @. We, 
therefore, consider a broad range of models. Since most 
the neutrinos are produced near the threshold, only the 
height of the EBL peak near 1 /im affects the results. 
At those wavelengths, the model of Ref. [22[ predicts a 
higher photon density than most other models. At the 
lower side of the range for 1 jzni EBL density are the 
models of Ref. [23[ and Ref. [25j. We show the spectra 
for these three models in Fig. [T] 

While AGNs are widely expected to be able to accel- 
erate cosmic rays to very high energies, little is known 
about the spectrum of cosmic rays produced by a typ- 
ical AGN. We have assumed the following form of the 
spectrum: 

j(E) (x E- a exp{-E/E max ) exp(-^ min / E). (2) 

The results do not depend on the lower energy cutoff, but 
the required source power does. We used E m i n = 10 13 eV, 
and we explored different values of -E m ax and a. The best 
fit to the IceCube flux (without overshooting the diffuse 
cosmic-ray and gamma-ray backgrounds) was obtained 
for a = 2.6, £ max = 3x 10 17 eV. We note that the 
corresponding gamma factor of a proton at the site of 
acceleration is close to the maximal value obtained in 
some detailed simulations 

The contribution of distant sources depends on their 
evolution with redshift. Following the analysis of 
Ref. |26j . we parameterize the source density evolution 
as 

(l±z) m , < z < zi 



Here m,zi,z%, and k are parameters obtained from fit- 
ting the observational data; they take different values for 
different AGN X-ray luminosities L x . Hasinger et al. j2(| 
obtain the following parameters based on observational 
data: 



Lx,erg/s 


10 42.5 


10 43.5 


10 44.5 


10 45.5 


m 


4.0 ±0.7 


3.4 ±0.5 


5.0 ±0.2 


7.1 ± 1.0 


Z\ 


0.7 


1.2 


1.7 


1.7 


Z2 


0.7 


1.2 


2.7 


2.7 


k 


-0.32 


-0.32 


-0.43 


-0.43 


Iff 1 n40 erg 
W P>W »MpcS 


7.0 


6.0 


1.3 


0.22 



TABLE I. Evolution parameters for AGN with different values 
of the X-ray power L x inferred from observational data [2(| 
are shown in the upper part of the table. The required power 
per unit volume W p of cosmic rays with energies E p > 10 13 eV 
was calculated under the assumption that an average AGN is 
described by one of these evolution models. 

We will consider all of these types of redshift evolution 
because one does not know whether the X-ray luminosity 
is well correlated with the power of cosmic ray emission. 

For each neutrino flavor we calculate the expected 
number of events in the energy interval of interest by con- 
volving their predicted spectrum with the experimental 
exposure given in Ref. [l[ . The overall flux normalization 
is chosen on the basis of the following criteria: (i) the pre- 
dicted average total number of neutrino events N v in the 
energy range 0.4 PeV < E < 6PeV must be as close as 
possible to the observed value N v = 2 (68% CL interval 
around 2 is shown in Fig. [TJ; (ii) the Poisson probabil- 
ity to observe at least 1 event above 6 PeV in the model 
must be less than 0.68, that is N% p < 1.14; (hi) diffuse 
photon flux should not exceed the Fermi upper bound; 
(iv) the predicted cosmic ray flux should not exceed the 
observed flux, for which we use the KASCADE-Grande 
results 
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Z\ < z < z 2 
z > z 2 



(3) 



The results of our numerical calculations are shown in 
Fig.[TJ As one can see from the figure, neutrinos produced 
in interactions of cosmic rays with background photons 
can account for the observed neutrino flux reported by 
IceCube collaboration in the case of strong evolution and 
high EBL of Ref. 0. 

The energy requirements per source are consistent with 
what is expected from AGN. For each of the models 
shown in Table Q] and in Fig. [TJ we calculated the emis- 
sivity at z — in cosmic rays with energies above -E m i n = 
10 13 eV. The results vary from 2 x 10 39 erg/s/Mpc 3 to 
7 x 10 40 er g/s/ Mpc 3 . Assuming the AGN density of 
10~ 5 /Mpc 3 [29|] , one obtains an individual AGN lumi- 
nosity of Lo ~ 10 44 erg/s for the lower end of the above 
range. This is a reasonable luminosity, which corresponds 
to the Eddington mass of 1O 6 M0. (AGN jets can ex- 
ceed the Eddington limit, but, in our case, the average 
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E, eV E, eV 

(c) (d) 

FIG. 1. Predicted spectra of PeV neutrinos (red lines) compared with the flux measured by the IceCube experiment 0]. 
The IceCube data points (red) are model-dependent 68% confidence level flux estimates obtained by convolving the IceCube 
exposure with the predicted neutrino spectrum. The predicted spectra are shown for the sum of three flavors; each flavor 
contributes, roughly, 1/3. The solid, dashed, and dash-dotted red lines correspond to the EBL models of Refs. [22l ]. [23| . and |25|. 
respectively. The evolution parameters for each plot are listed in table Q] for (a) L x = 10 42 5 erg/s, (b) L x = 10 43 5 erg/s, (c) 
L x = 10 44 ' 5 erg/s, (d) L x = 10 45 5 erg/s. In all cases, we assumed the proton spectral index a = 2.6 and the maximal proton 
energy -E max = 3 x 10 17 eV. Also shown are the predicted gamma ray (lower curves below 10 TeV) and cosmic ray (upper curve) 
fluxes. The cosmic ray data points above 10 PeV are based on KASCADE-Grande the diffuse gamma-ray background 
data points below 1 TeV are due to Fermi [2^| . 



AGN is expected to operate well below the Eddington 
luminosity.) This is also consistent with the analyses of 
Refs. 

Future results from IceCube may help constrain mod- 
els of cosmic ray acceleration in AGN. We note that 
cosmic ray flux provides a stronger constraint than the 
diffuse gamma-ray background. Composition measure- 
ments based on the data of KASCADE-Grande 27| are 
subject to large uncertainties in the Monte Carlo simu- 
lations, especially in the energy range of interest to us. 
Furthermore, local galactic magnetic fields can affect the 
flux and composition of cosmic rays with energies be- 
low 10 17 eV (and even those with higher energies 30]), 
making it difficult to connect the locally measured com- 
position to that of extragalactic sources. Therefore, we 



used the total cosmic ray flux as the upper bound. 

In summary, we have examined the recent observations 
of the IceCube experiment in light of the model that ex- 
plains the spectra of distant blazars by secondary gamma 
rays produced in cosmic-ray interactions along the line of 
sight [4-ll2|. The same interactions result in a neutrino 
spectrum peaked at 1 PeV We have shown that dis- 
tant AGNs can generate a spectrum of neutrinos with a 
peak at 1 PeV and with a flux that is consistent with the 
IceCube results. 
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